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Abstract 

In this paper, employing the modified dispersion relation, we have derived the general mod¬ 
ified Friedmann equations and the corresponding modified entropy relations for the Friedmann- 
Robertson-Walker (FRW) Universe. In this setup, we find that when the big bounce happens, 
its energy scale and its corresponding modified entropy behavior are sensitive to the value of k. 
In contrast to the previous work with k = 0, our work mainly demonstrates that the bouncing 
behavior for the closed Universe with k = 1 appears at the normal energy limit of the modified 
dispersion relation introduced, and when bouncing phenomenon is in presence, its modified entropy 
is just equal to zero. Surprisingly, when k = —1, the bouncing behavior is in absence. 
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I. INTRODUCTION 


General relativity plays an important role in deeply linking the geometric structure of 
spacetime with the distribution of matter sources. This remarkable theory provides a pow¬ 
erful perspective to understand the nature of gravity and describe the evolution of the 
Universe. With the cosmological principle, the standard cosmological model can be derived 
successfully from general relativity. However, the cosmological singularity problem cannot 
be avoided effectively in this model. To solve this singularity problem, theorists widely turn 
to appealing for the quantum gravity. Unfortunately, up to now, one has not yet established 
a complete and consistent quantum theory of gravity, which results in that the singular¬ 
ity problem cannot be solved at hrst principle level. As a result, people need to seek for 
some important quantum properties in string theory and loop quantum gravity as important 
mechanisms to relieve or solve this problem at a phenomenological level. 

Nowadays people generally believe that the combination of general relativity and quan¬ 
tum mechanics will provide a fundamental minimal length or the maximal energy, which 
is referred as a feature of quantum gravity. The developments of string theory and loop 
quantum gravity have greatly intensihed this belief. In particular, based on Loop Quan¬ 
tum Gravity there have been great interests in some modihcation of dispersion relations 
I-3|, and Generalized Uncertainty Principles(GUP)y, have been taken into account in 
the literature on string theory [t] and on noncommutative geometry 8]. The forms of 
the modihed dispersion relation (MDR) and ones of GUP can be viewed as the alternative 
mechanisms to handle the black hole problems, such as the singularity, the “thermodynam¬ 
ics”, the evaporation and the “information paradox” and so on. In the recent years, the 


implication and application of the effects of quan 


attentions 


32l | ( for recent reviews we refer to 


um gravity have attracted a great deal of 
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35| 


In the semi-classical limit an intuitive scenario was proposed to avoid the cosmological 
singularity problem in standard cosmology by a big bounce 36|, l37| that was implemented 
by modifying the standard Friedmann equation. Recently, this scenario has been extensively 
investigated in various approaches 38l-l48| . In particul^ based on the assumption that the 
gravitational theory can be described by entropy force 49|], which has been widely interested 
by many theorists [40|, l50l-l63|. the authors in proposed that an interesting setup was 


that the modihed dispersion relation prevented the energy density of the matter contents 
from diverging at high energy level such that the modihed Friedmann equation with the 
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bouncing effect can be derived successfully from the Clausius relation, where the singularity 
of the corresponding spacetime was free. 

In the previous work the authors had discussed the relevant bouncing behavior 
and the corresponding entropy behavior for spatial flat universe by the modified dispersion 
relation. Now we wonder whether similar interesting behaviors for the cases with non¬ 
vanishing k will happen. We expect to provide an answer to this problem in a parallel way. 
It turns out that using the modihed dispersion relation in 40| to prevent the energy density 
of the matter from diverging, we can obtain the general modihed Friedmann equations 
and the corresponding modihed entropy relations for the FRW Universe. When A; = 0, 
our results naturally reduce to that considered in 40|. In particular, we hnd that the 
modihed Friedmann equation contains bouncing behavior and the corresponding entropy 
is non-negative for the closed universe. More precisely, its bounce appears at the normal 
energy bound of the MDR, namely E = This implies that when bouncing phenomenon 
is in presence, its entropy is just equal to zero, which is quite diherent from the result in the 
previous work In addition, we have noted that when k = —1, the bouncing behavior 
is in absence. Our results have strongly suggested that k plays an important role on the 
bouncing behavior of the Universe. 

Our paper is organized as follows: in Sect. II, we briehy review on the description of 
entropy force for the Friedmann-Robertson-Walker universe; in Sect. Ill, our task is to deter¬ 
mine the modified universe equations with the MDR in detail and discuss the corresponding 
entropy relations; in Sect. IV, we discuss the bouncing effects of the universe with various 
values of k in detail; Sect. V gives conclusions and discussions. In Appendix, we present a 
detailed calculation to give out the reason why the bouncing behavior for the closed universe 
occurs at = rjlp. Here we have set h = c = 1 for convenience. 


II. THE ENTROPY FORCE DESCRIPTION 

In the section, we will briefly recall some ingredients for deriving the Friedmann- 
Robertson-Walker universe from the relation of the entropy force. Let us begin with in¬ 
troducing the Friedmann-Robertson-Walker metric which is 

ds'^ = -dt^ + 2 + ( 1 ) 

1 — kr-^ 
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where k values can be taken —1,0 and 1 which correspond to the open, flat and closed 
Universe, respectively, and a{t) is the scale factor of the Universe. This spacetime geometry 
implies that the apparent horizon 6^ can be defined as 

1 


rA = 




k 


( 2 ) 


In analogy with black hole thermodynamics, for the apparent horizon we can correspondingly 
define its entropy as 

An , , 

( 3 ) 




as well as the temperature 


4G G{m + 


T = 


2nrA' 

where A = 4:nr\ is the apparent horizon area. In the viewpoint of entropy force, utilizing 
the Clausius relation, SQ = T6S , and integrating the both sides of this relation, one can 
naturally obtain Friedmann equation 


k 8nG 


H^ + ^ = 


-Pu 


( 5 ) 
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where pt is total energy density. Here we have successfully derived Friedmann equation from 
the Clausius relation via assuming that our universe is homogeneous and isotropic, and 
treating one as a thermodynamical system in equilibrium. It is obvious that the singularity 
still exists. Thus, in the next sections, after taking a specific MDR into account, we expect 
to apply the corresponding entropy force prescription to construct the modified universe 
equations with any spatial curvature, which may include the bouncing universe behavior, so 
that the singular problem can be free. 


III. MODIFIED UNIVERSE EQUATIONS WITH MODIFIED DISPERSION RE¬ 
LATION 


In this section, we mainly concentrate on the modified dispersion relation how it causes 
the modification of the FRW Universe equations with a spatial curvature k in semi-classical 
limit. Here we will introduce the quantum gravitational effect at the phenomenological level. 
Let us start with the modified dispersion relation. Based on Doubly Special Relativity, a 


general modified dispersion relation proposed in 


11 


12| has a following form. 


E^f\l,E)-p^g\l,E) = m\ 


( 6 ) 
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where E, p and m are the energy, momentnm and rest mass of the particle, respectively, 
= SnG. The above MDR has been widely used to study black hole physics and rainbow 


universe 


15-19 


3,I3- 


At low energy limit, namely IpE <C 1, one naturally requires that 
the two general functions f{lpE), g{lpE) ~ 1, so that the above dispersion relation reduces 
to the standard Einstein energy-momentum relation. The authors of the previous work 40| . 
in fact, have chosen the two specific functions for the modified dispersion relation 

sin^ {r]lpE) 


then Eq.( 


nipE) = 

gives rise to the following form. 


g (IpE) = 1, 


(7) 


^ sm{r]lpE) = (8) 

where g is dimensionless parameter which is bigger than zero. They have omitted the 
negative energy branch. The right side of the equation is not less than zero, so that the 
value of energy is in the range [0,^]. Obviously, this MDR can be separated into two 
branches, namely [0, and ^], respectively. It is easy to hnd that in the different 
branches, the monotonous behavior of their momentum is remarkably distinguishing, as the 
energy is increasing. More explicitly, for the particle energy in branch [0, its momentum 
p is increasing as the energy is growing, while for the particle energy in branch [^j-, ^], 
p has a strangely and monotonously decaying behavior. Thus we call the latter branch as 
the anomalous energy region. In particular, for the rest mass m = 0, in between the two 
branches, when particle energy passes through E = its momentum reaches a peak 
value, namely^; whereas when E = ^, its momentum p vanishes. In order to introduce 
the bouncing effect in the Universe, we need to establish further the variation relation of the 
above equation ([H]) that has a following form 


5E = ± ^ ^ =6p, (9) 

\/l- sin^{glpE) ^{p^ m?) 

where the plus sign “ -|- ” and the minus sign “ — ” stand for choosing the energy value in 
the range [0, and respectively. For the high energy limit or massless particle, 

the above equation can be expressed as 


5E 


± 


\/l - {glppf 


6p. 


( 10 ) 


Later, we will see that the above relation plays an important role in deriving the modified 
universe equation which may relieve the singularity problem. When the modified dispersion 


5 























relation is taken into account, the modified Friedmann equation can be derived via an 
integration relation in 62 | 


StiG 


P=-^ I S\A){^fdA, 


AG. 


G 


A 


( 11 ) 


where 5"(y4) = ^ is related to the modified dispersion relation, while for ordinary case 
it is a well-known constant, namely Now let us consider how such modified ratio of 
the entropy to the area on the apparent cosmological horizon is derived out, and then how 
it governs the modified Friedmann equation such that the singularity of the universe can 
be free. We treat our universe as a thermodynamical system in equilibrium and suppose 
a quantum process in which a single massless particle with energy E escapes thr oug h its 
apparent horizon. This causes correspondingly the minimal entropy change Q m Q, 


which should satisfy Clausius relation. 


5Q 


dSmin = dS = — = ± 


271 


4:Tr{vlp)^ 

A 


( 12 ) 


here we have used the identification relations 60 


Heisenberg uncertainty relation 6p 


1 1 
rs-/ —— r\j - 

OX VA 


20 |. 


E ~ 6E as well as p ~ 5p, and the 
Moreover, we have also admitted a 


physical fact that for the particle in the quantum process there is an intrinsic uncertainty 
position as its Compton wavelength that is identified with the apparent horizon of the 
universe. In this quantum process, correspondingly, one assumes that causing the minimal 
area change of the apparent horizon is dA^in = = 8"^^ j40|(and references therein). 

Identifying dA with 6Amin and using Equation ([12]), we have 

1 1 


dS _ dSmin 
^ dA~ 6 Amin 


AG 


A 


(13) 


Let us first consider the plus case. Putting this equation into the integration relation (ITT]) , 
we can obtain the modified Friedmann equation. 


8tvG 


-P = 




AiT'qH'l 


A 


+ Cl, 


(14) 


when <C A, Eq. flTTl) should reduce to the standard Universe equation, which is 


thus the integral constant can be figured out 


Ci = 




A 


(15) 
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So the modified kinematic Friedmann equation can be determined as 


SttG SttG a 


-Pt = 


3^3 7]HI 


(1 


471-772/2 

^ 4 “' 


( 16 ) 


Making both sides of Eq. ffT 6 |) square and simplifying it, we can obtain an alternative form 
as 

h n, 

(17) 


2 , k SttG Pt 

H + — = 1-, 

a2 3 Pc 


where pc = ^4. With the use of the continuity equation pt + 3H{pt + Pt) = 0, the modified 
dynamical Friedmann equation can be calculated out 


- = ^Pi(l - -)-i^G(p, + P,)(l - 2^), 

(X o pc pc 

where Pt is the total pressure. 

The modified entropy-area relation for energy range [0, 1 ^] can be expressed as 


(18) 


A I 47rr7^ 77772/2 / / A 

4gV^ a ^ G 477772/2 + Y 47r772/2 


1] + Di, 


(19) 


when 77 —>• 0 , the modified entropy-area relation reduces to the normal result, namely 
which can give rise to the integral constant Di = 0 . 

Similarly, for the energy range [ 2 ^) 7 ^] we can work out the modified Universe equation 


SttG 


-P = 


3 " 772/2 V A 


1 - 


477772/2 


+ G 2 . 


( 20 ) 


Here we require that when A = 477772/2 corresponding to the energy value Eq. fl20|) can 
match smoothly with Eq. flTT)) so that G 2 should be equal to Gi, namely G 2 = Ci = -4j7 


A 

3 • 


In the parallel way, with the use of the above equation and continuity equation, it is not 
hard to reproduce the modified universe equations which are the same as Eqs. fll7l) and flTSD . 
respectively. Similarly, the corrected entropy relation also becomes 


A 1. 477772/2 


Sm — ——\l 1 — 


4G 


77772 /, 


2/2 

'P 


A 


G 


In). 


A 


47r7]^P 


+ 


A 


477772/2 


1] -|- D 2 . 


( 21 ) 


Here we also require that the modified entropy-area relations in both Eqs. flT^ and fl2T]) are 
the same, when A = 477772 / 2 . Thus we can obtain the integral constants satisfy a relation 
D 2 = Di = 0. It is easy to check that when A = the value of entropy-area vanishes, 

namely Sm = 0 . 
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IV. THE BOUNCING BEHAVIORS IN THE MODIFIED UNIVERSE EQUA¬ 
TIONS 


In this section, we mainly explore the bouncing behavior for the modihed Friedmann’s 
equations. In the bouncing universe scenarios, the most appealing feature is no initial 
singularity of the universe. In these bouncing paradigms, ones usually argue that there 
should be a nonsingular connection between the two distinguishing phases that are the 
contraction phase and the expansion phase in order. More explicitly, the evolution of the 
cosmological scale factor from the contraction phase with a < 0 to the expansion phase 
with a > 0, goes through its minimal value with non-vanishing value referred as a critical 
point. In such process, some effects (quantum effects) that become the dominant ingredients 
would prevent the universe from collapsing into a singularity and then drive our universe 


ogical scale reaches 


40[. Thus, we can 


to accelerate expansion. As a result, at the critical point the cosmo 
the minimal nonzero value (oq > 0) which has Ho = 0 and ho > 0 
reasonably view them as the bouncing conditions. After taking the hrst bouncing condition 
into account, namely Hq = 0, solving reversely modihed Friedmann equation ffT7|) gives rise 
to 


1 / 

= ( 22 ) 

which should meet the constraints pt± > 0 and ta > 0. Here we have denoted the corre¬ 
sponding scale factor as oq when the Hubble parameter vanishes. Making use of Eqs. flTT)) 
and (12^ . the modihed dynamical Friedmann equation is correspondingly represented by 


Qq 

Oo 


X±A7rG 

% 



kr]‘^lp 

al 


{Pt± + Pt), 


(23) 


where the signs “ ± ” in the above equation correspond to pt±, respectively. There are two 
points what we would like to emphasize. Firstly, pt- is an ordinary energy density in the 
energy scale region [0, 2 ^]) while pt+ is anomalous energy density in [ 2 ^, ^]- When taking 
energy density relation in fl2^ to calculate the cosmic acceleration in between the contracting 
phase and the expansion phase, we have to use the corresponding relation in fl2^ . Secondly, 
when k ^ 0, this means there is an undetermined degree of the freedom, namely the scale 
factor oq. Now we are going to divide them into three cases, namely k = —1, A: = 0 and 
k = 1, respectively, and check whether they really satisfy the constraints and the other 
bouncing condition or not. 










• For k = —1, when requiring this case to satisfy the first bouncing condition, one 
easily finds that becomes an imaginary number, which is not consistent with the 
constraint > 0. Thus for this case the bouncing solution does not exist! 

• For k = 0, we have 

Pt+ = pc;pt- = 0. (24) 


When our Universe reaches high energy bouncing critical point which corresponds to 
the total energy density is rather pt+ = Pc than pt- = 0, this gives out the dynamical 
Universe equation which is 


ho 

Oo 


— AttG^Pc + Pt) > 0 , 


it is iust the result 




(25) 


• For k = 1 the total energy density equation naturally reduces to the following 
form 

It seems that there are two bouncing solutions. Based on the definition of the apparent 
horizon and the MDR, we can hnd the scale factor is no less than rjlp. In Appendix, 
according to the bounce conditions, we can show that the bouncing effect of the 
Universe is presented at the minimum scale value that means the momentum reaches 
the maximum value corresponding to the critical point E = such that in fact 
the above equation (l26ll can only take the same value, namely pt- = pt+ = y with 
Oq = After the big bounce of our Universe happening, the total density will 

decrease monotonously with our Universe expanding. When the total energy density 
climbs up the value pt- = pt+ = y, Eq. fl2^ can be rewritten as 


Oq 


ao 


> 0 . 


(27) 


Thus for this case the bounce of the Universe can exist in the normal energy region 
[0, 2 ^] without anomalous negative entropy! 


Here it is worth noting that for the case /c = 0 the bouncing solution occurs at the crit¬ 
ical point with maximum energy value but with negative entropy in For k = 1 
the bouncing phenomenon can be presented at the energy value corresponding to the 
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maximum momentum ^ . This means that the modihed entropy-area relation is rather 
Eq. flT^ than Eq. flMll . More precisely, when the universe performs its bouncing behavior in 
high energy limit, its corresponding entropy-area relation is just equal to zero. After the big 
bounce, its modihed entropy will increase monotonously with our Universe expanding. For 
k = —1 the bouncing one is absent, since it can not meet the constraint condition > 0. 


V. CONCLUSIONS AND DISCUSSIONS 


In this paper, along with the spirit in ^ , assuming that both the Clausius relation and 
the modihed dispersion relation keep ehective on the apparent horizon, we have derived 
the general modihed universe equations and the corresponding modihed entropy relations 
for the FRW Universe. From the modihed universe equations, we hnd that the bouncing 
behavior is presence for the closed Universe, besides the spatial hat Universe, so that the 
cosmological singularity problem can be solved for these cases at the phenomenological level. 
In particular, in contrast to the bouncing solution of the spatial hat universe with A: = 0 
to occur at the maximum energy value ^ of the MDR ([8]), but with negative entropy in 
j^ . we have shown that the bouncing behavior of the closed Universe with A: = 1 explicitly 
appears at the energy value without negative entropy. More precisely, when the closed 


universe performs the bouncing behavior at the critical energy point the corresponding 
entropy-area relation is just equal to zero. Surprisingly the bouncing solution of A: = —1 is 
absent, since it can not satisfy the constraint condition > 0. 

However, it is well-known that the entropy relation of the Universe is an open question in 
the bouncing or cyclic universe scenarios, if we respect the second law of the thermodynamics 
saying that the total entropy of the universe never decreases as time flows. When the closed 
Universe evolves from the contraction phase to the expansion phase, Eq. flT^ suggests that 
the corresponding entropy-area relation has diminishing behavior and increasing behavior, 
respectively. These facts tell us that the entropy behavior in the contraction phase does not 
agree with the second law, but one in the expansion phase agrees with the second law. Thus, 
in analogy with the previous bouncing or cyclic universe scenarios, the well-known entropy 
problem that violates the second law is still in presence in our scenario. A deep understanding 
on how to link the diminishing behavior of the entropy in the contracting phase and the 
second law of the thermodynamics is still lacking. There might be an approach to solve 
this problem via considering the generalized second law of thermodynamics which has been 
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discussed in many universe models 64l-l67l|. Under complying with the second law of the 


thermodynamics, it should be interesting to seek for some new modihed dispersion relations 
and new physical mechanism to overcome or solve the singularity problem of the Universe. 
This aspect is left for future works. 
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Appendix 


In this Appendix, we present a detailed calculation which demonstrates that the bouncing 
solution for the closed universe can appear at the minimum of scale factor. To obtain the 
scale factor size at the critical bouncing point, we need to take into account the behavior 
of the total energy density. When our universe evolves to around the critical point in high 
energy limit in the contracting phase, we are interested to consider the behavior of total 
energy density with the constant state parameter u that takes the values in the range 
(—|,oo). Since in this case the cosmological constant density is very small and can be 
totally ignored, the behavior of the universe is governed by the dominated-content with the 
constant state parameter u. Thus for any given the state parameter value u in the allowed 
region, we can identify the dominated-content density with the total energy density. From 
the continuity equation, we can get the following density relation 


C 

a3(i+^) ’ 


(28) 


where C is the integral constant. From the Friedmann equation (lT7|) . it gives rise to 

(§)" = ^aV<(l - - 1 > 0. (29) 

at 3 Pc 

When the above equation is equal to 0, it means the Hubble parameter H is zero which 
is just the hrst bouncing condition. Substituting Eq. fl28l) into the above equation, we can 
obtain 

36^4+60. _ i2/2c'a3+3<^ + < 0, (30) 
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for pt < Pc and a > rjlp this equation should be valid. Now let us determine the integral 
constant C. When putting a = rjlp into Eq. (l30|) . we can easily get 

rjHliC - < 0, (31) 

on the other hand, the above equation itself has 

> 0, (32) 

it is worth noting that this is the constraint exactly what we are looking for. This rela¬ 
tion plays an important role in deciding the critical point. As a consequence, combining 
Eqs. (lM]) and (l3^ . we have 

C = (33) 

Plugging it into Eq. fl30|) . we hnd 

^4+6a; _ < 0. (34) 

Eqs. flTT]) and fl5^ suggest that Oq = plp is a root of the Friedmann equation flTTD satisfying 
the hrst bouncing condition. It is easy to check that this root saturates totally the other 
bouncing condition, as showed in Eq. fl27)) . However this story is not yet end so far, since 
Eq. fl3T)) generally has other roots. Therefore, in the next step, we have to show that other 
roots can not satisfy the second bouncing condition, otherwise the bouncing behavior occurs 
at the other roots. Here we provide a general trick to discuss this problem without specihc 
state parameter u. For convenience we dehne a function 

/(a) = (35) 


In order to obtain the monotonicity of this function, we need to work out the derivative of 
the function with respective to the scale factor a, which is 

^ ^ = (4 + 6u,)a^+=>“[a>+^ - = 0, (36) 


for o; > — I and a > 0, we can obtain 


a = 


3 -|- 3a; 
*-2 + 30 ; 


I 1+3“ pi 


P- 


(37) 


Thus, we can easily check that /(a) decreases monotonously, when the scale factor a in 
the range [plp,a]; f{a) increases monotonously, when the scale factor a is outside a. Since 
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«-o = vh is a root of /(a) = 0 and /(a) < 0 , if it has other roots, then they must be greater 
than a. 

Now we turn to considering the role of the second bouncing condition. In order to seek 
for the bouncing solution in normal energy region, we should consider the minus cases in 
both Eqs. fl22|) and fl2^ . Substituting pt- = \pc0- ~ \J^ ~ Eq.([23]), we find 


iio _ 3(1 + 0 ;) 2 + 3a; 3(1 + cj) [ 

Requiring that the above equation is greater than zero, we can derive an important constraint 

2 (3a; + 2)Wll 

° (3a; + 3)(3a; +1) 

It tells us that if any one root Uq meets the above requirement, then the bouncing behavior 
can appear at the root. It is easy to check that Oq = plp satishes this constraint. In following 
we will prove that the other roots do not satisfy the requirement. For the purpose, we can 
define a function which has a following form 



/ \ (3a;+3) (3i^+l) 

9{uJ) = ^ 

. . 4+6u; 

(3a; + 2 ) 3^+1 

^ . 3+3a; . ” ? 

( 3 a; + 3 ) si^+i ( 3 a; + 1 ) 


(40) 


then taking it logarithm, we have 


y{iS) = (3a; + 1) In 51 ( 0 ;) 

= (3a; + 3) [ln(3a; + 2) - ln(3a; + 3)] + (3a; + 1) [ln(3a; + 2) - ln(3a; + 1)]. (41) 


To determine the monotonous behavior of this function, we can find its derivative with 
respective to the state parameter u. 


dy{u}) 

doj 


3 In 


(3a;+ 2)2 
(3 + 3a;)(l + 3a;) 


> 0 . 


(42) 


Thus the function y{uj) increases monotonously in the range (—|, cxa). Now we work out the 
behavior of the function at the infinity. Making Eq. flTT]) take limit, we can have 


lim y{uj) 


Q I O 

lim [(3a; + 3) (In —-) + (3a; + 1) (In 

o;—>00 oCJ o 

,• r 1 1 

lim-h ... 

oj^oo 6 + 60; 2 + 60; 


3a; T 2,, 


0. 


(43) 
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So, the function y{u) is less than zero in (—i,oo), which means that the function g{u) is 
less than 1. This tells us that although the other roots can exist, they can not saturate the 
second bouncing condition, namely do > 0. Thus, the bouncing solution appears at uq = rjlp. 
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